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Abstract 

We present the full one-loop radiative corrections to pair production of neu- 
tralinos in e + e~ collisions within the Minimal Supersymmetric Standard Model. 
Particular attention is paid to the definition of weak and QED corrections. The 
non-universal QED corrections are extracted by subtracting the initial state radia- 
tion. We give numerical results for two different SUSY scenarios for e + e~ — ► X1X2 
and e + e~ — ► ^2X2- The weak and QED corrections are up to several percent or 
even higher and need to be taken into account at future linear collider experiments. 



1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) [Ij, one has two charginos xt 
and xti which are the fermion mass eigenstates of the supersymmetric partners of the 
W ± and the charged Higgs states H 1 2 . Likewise, there are four neutralinos Xi~X% which 
are the fermion mass eigenstates of the supersymmetric partners of the photon, the Z° 
boson, and the neutral Higgs bosons 2 . Their mass matrix depends on the parameters 
M, M', /i, and tan/5, where M and M' are the SU(2) and U(l) gauge mass parameter, 
and tan j3 — ^ with t>i >2 the vacuum expectation values of the two neutral Higgs doublet 
fields. If supersymmetry is realized in nature, charginos and neutralinos should be found 
in the next generation of high energy experiments at Tevatron, LHC and a future e + e~ 
collider, . Especially at a linear e + e~ collider, it will be possible to perform measure- 
ments with high precision j2j 0] . In particular, it has been shown in j2] that the masses of 
charginos and neutralinos can be measured within an accuracy of Am^±,o =0.1 — 1 GeV. 
It is therefore obvious that such a high precision requires equally accurate theoretical 
predictions. Despite the complexity, for some SUSY processes the full one-loop correc- 
tions have already been calculated: for e + e~ — > XiXfi hj = 1)2, in p[|, for e + e~ — > klj, 

I = e, fi, i,j = L, R in [Sj, e + e~ — > fa fj, f = q,l,u (including the third generation) in 
[HllZj. As to decays, the full one-loop corrections were calculated for q { — > qx®, 2 = 1 — 4, 

and fa — > qXk) ^ = 1, 2, in jH], and for the decays A — > fxf 2 , H — * /i^4o m [0], where A 
is the pseudoscalar Higgs particle. All these calculations have shown that the corrections 
are important for precise predictions of cross sections, branching ratios and asymmetries. 
In this paper, we present the calculation of the complete one-loop corrections to the neu- 
tralino production e + e~ XiXp ^, j = 1 — 4. 

For the calculation of higher order corrections, renormalization of the MSSM is neces- 
sary. For this purpose, one has to employ appropriate renormalization conditions, or 
equivalently, one has to fix the counter terms for the SUSY parameters. In this paper, 
we adopt the on-shell scheme for the chargino and neutralino system of [TO] . Equivalent 
methods were developed in [TT) |H] . The schemes only differ in the fixing of the counter 
terms of the parameters M, M' and \i. Hence the meaning of these parameters is differ- 
ent at loop-level. The schemes, however, yield the same results for observables as masses, 
cross sections, widths, etc. . 

Starting from the tree-level in section |2l we outline the calculation of the one- loop correc- 
tions in section El discussing the renormalization both of the SUSY and SM parameters. 
The process-independent corrections to the neutralino mass matrix are included in an 
improved tree-level. Particular attention is paid to a proper definition of the weak and 
QED corrections as they latter play an important role. In section we represent a de- 
tailed numerical analysis for e + e~ — > X1X2 an d e + e~ — > X2X2 f° r a higgsino and a gaugino 
scenario for Xi an d X 2 - Conclusions are given in section |3J 
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2 Tree-level 



In the MSSM the neutralino sector is specified by the gaugino mass parameters M and 
M', the higgsino mass parameter /i and the Higgs mixing angle tan/3, all appearing in 
the neutralino mass matrix (in the bino, jy 3 -ino, if 12 -ino basis) 



Y 



/ M' — mz sin 6 W cos ft m z sin 6 W sin ft \ 

M mz cos 6*iy cos /? — mz cos 0^ sin ft 

—mz sin 0^ cos ft mz cos 0yp cos ft — /i 

V sin W sin /3 — m z cos 6^ sm /5 — A* / 



With the unitary matrix iV, which diagonalizes the mass matrix Y 

diag(m^o, m^o, m^o, m%o) = iV*FiV^ , (2) 

we can rotate from the gauge eigenstates ip® = —i\ 3 ,ip H - 1 ,ip H - 2 )j to the neutralino 

mass eigenstate basis x° = Nijip®. 

At tree-level and neglecting the electron mass in all Yukawa couplings the production 
process 

e + e~ ^x-X°j (hJ = 1,2,3,4) 

contains contributions from the Feynman diagrams shown in Fig. ^ The direct s-channel 
due to the Z° exchange and the crossed t- and w-channel due to the et,R exchanges. 




Figure 1: Tree-Level 

From the interaction Lagrangian 

£z°ee = g —Zle^[C L P L + C R P R ]e, (3) 

COS &w 

= gtfeP R x°e L + gf t R eP L x°e R + h.c., (5) 
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we obtain the couplings 



Cl,r 


= J3L,R _ 




= -0" R 






ft 






2 



0w, I* L = -\, ^ = 0, (6) 
±N l3 N; 3 + ±N tA N* A , (7) 



(8) 



3 One-loop corrections 

The radiative corrections to the neutralino pair production include the following generic 
structure of one-loop Feynman diagrams: The virtual vertex corrections Fig. El the cor- 
rections to the tL,R and Z° propagators Fig. El and the box graph contributions Fig. |H 
The notation F, V, and S stand for all possible fermion, vector and scalar particles in 
the MSSM, respectively. U denotes the FP ghosts. Diagrams with loops on the external 
fermion lines are included in the definition of the counter terms as wave function correc- 
tions. In this work, the complete set of Feynman graphs is calculated with help of the 
packages FeynArts and FormCalc [12 . We implemented our renormalization procedure 
into these packages. For a proper treatment of the appearing UV divergencies, counter 
terms are introduced in the on-shell renormalization scheme. To preserve supersymme- 
try, the used regularization scheme is dimensional reduction (DR). The loop graphs with 
virtual photon exchange also introduce IR singularities. Therefore, real photon emission 
has to be included to obtain a finite result. 

a ™rr {e+e - _> ^0-0) = a ren {e+e - ^ ~0 ~0) + ^ & - ^ ~0 -0^ (g) 

For the numerical analysis, we have also used the programs Loop Tools and FF [T5] . 



3.1 Renormalization 

Wave function counter terms 

In the prescription of the used on-shell renormalization scheme all involved fields get the 
following shifts to obtain the so-called wave function corrections. 

x° (% + m,p L + iszSKifS , ( A ) _ ( i + isz L i + (10) 
z. . (1 + i< z„)z. + ^, (f^( 1 + l SZ[ 1 + l s4 )(l)- (») 

with the definition of the renormalization constants 
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e+ 



e+ 




*- 



+ crossed graphs 




+ crossed graphs 



s/v x° 

+ crossed graphs 



T e. 



5/V 




x° 



+ crossed graphs 



Figure 2: Generic Vertex Corrections 





+ crossed graphs + crossed graphs + crossed graphs 



Figure 3: Generic Propagator Corrections 



SZ ZZ = -»n Z z(m|), 5Z Zl = 2 ^ n f (0) , (12) 



*z£ = -KIlLKj, 6Zi = -KriLKj, (13) 
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crossed graphs 



crossed graphs 



Figure 4: Generic Box Corrections 



sz 1 



5Z 



6Z. 



5Z 



R 



3? 



- n L (m 2 f ) - m 2 (fl L (m 2 f ) + II R (m 2 )) 



2m 



-{Tl SL (m}) 



IT SR (m 2 )) 



-m f (U bL (mj) + IP K (m 2 )) 



(14) 



- n-(mio) - mio(n^(m|o) + n£(m£o)) + — 



:nf (m| ) - nfK„)) 



-m,o(nf(m|o)+nfKo)) 



(15) 



m~oII^(m~o) + m^om^oUfAmio) + myonf^fm-o) + myoIIf/Vm-o 



rSR/ 



JZ^L <- i?) , = - J2) , 5Z£ = SZh(L <-> R) 



(16) 
(17) 



where n fe) (fc 2 ) = ^nf w (* a ) + *JWif y) (* 2 ) + i^nf^ 2 ) + P«ngJ5(P), il(m 2 ) 

^|2-n(A; 2 ) fc2 2 and Cjj = 2/(m| — wijq)- For the neutralinos it holds IT^(p 2 ) 



H^ L (p 2 ), because of their Majorana nature. Since we neglect 



Xi 



n^(p 2 ),nf /L ( P 2 ) 

the selectron mixing, no sfermion mixing angle need to be renormalized. 



Neutralino and sfermion mass matrix renormalization 

In the MSSM, the four neutralino masses depend on the SUSY parameters M', M, /i, 
and tan/? and the SM parameters mz and sin 9w As M, fi and tan/? also enter the 
chargino mass matrix, the renormalization of the neutralino, chargino and SM sectors is 
interrelated. Therefore, it is necessary to take into account radiative corrections to the 
X° masses and the rotation matrix. For the on-shell renormalization two different ap- 
proaches are essentially known in the literature, ^U] and [TT| |H] . Although the corrections 
in the neutralino masses are in general small, these shifts can lead to large effects near 
the threshold. It would be possible to adopt a renormalization scheme for each channel 
in such a way, that the two produced neutralinos are input parameters and do not obtain 
mass corrections. A threshold shift would thus be avoided, but this leaves us with the 
problem that the renormalized processes have unequal counter terms for different produc- 
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tion channels, which would lead to different meanings of the neutralino and chargino mass 
parameters. Here we use the on-shell scheme described in ^H]- We define an improved 
tree-level, where the pro cess- independent mass matrix renormalization is already included 
and separated from the residual weak corrections. Absorbing the finite correction AiVjj- 
to the rotation matrix Ny in the improved tree-level is equivalent to defining an effective 
coupling matrix + A . 

This yields the following counter terms for the neutralino mass matrix SYij and the rota- 
tion matrix SNij. 



SYaa 



5N i:i 



N H N nj Re [m^n^mjo) + m^<(m|o) + l&{m%) + nf n L (m|o ; 



Ln=l 



- £ (b~Zl - 6Z&) N kj . 
4 k=i 



(18) 
(19) 



The same renormalization prescription can be applied to the sfermion sector. Counter 
terms for the SUSY breaking masses Mq £ and Mfj ^ both entering the sfermion mass 
matrices, are introduced. Fixing Mq i in the down-type mass matrices results in a cor- 
rection to the up-type masses and mixing angles jHj. Hence, in our case, we have no 
additional corrections to the selectron masses. The correction to the electron sneutrino 
mass, which only appears in loop graphs, is of higher order and do not need to be con- 
sidered. 



Renormalization of the SM parameters 

Since we use as input parameter for a the MS value at the Z pole, a 
e 2 /(47r), we get the counter term P ITS] 



mm 



zJ\MS 



5e 1 e 2 

e (4?r) 2 6~ 



4E^H(A + iog§)+EE^(A + lo F <r 



■ f 

2 f o 2 
+4£(A + logJr 

k=l V " L 



■I i 



m=l 



171: 



E(A + log4^)-22(A + log^ 



(20) 



with Xf = m z V rrif < m z and x t = m t . Nq is the colour factor, N[, = 1, 3 for (s)leptons 
and (s)quarks, respectively. A denotes the UV divergence factor, A = 2/e — 7 + log47r. 
The masses of the Z boson and the W boson are fixed as the physical (pole) masses, 



6m 2 z = ^U zz {m 2 z ) 



bm 2 w 



Kfl^(m^) 



(21) 



and sin 2 6w is fixed by cos 6w = mw/mz- 
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3.2 Definition of weak and QED corrections 



As mentioned before, the full one-loop corrections become IR convergent if also real photon 
emission is included in the calculation. Because of these large additional corrections, it 
is desirable to treat the weak and QED parts separately. The easiest way to define pure 
"weak corrections" would be to separate off all Feynman graphs with an additional photon 
attached to the tree-level diagrams. However, in our case this cannot be done in a gauge 
invariant and UV finite way due to the selectron exchange channels. Another possibility 
would be to use the soft photon approximation [15], where only "soft" photons up to a 
maximal energy AE are included: <j weak = u^f 1 and a® ED = a . The weakness of this 
definition is the large AE dependence of the weak and QED components oc log The 
sum of both is, however, cutoff independent. Therefore, we extract the AE terms and 
the leading logarithms ^L e = - log rf? , caused by collinear soft photon emission, from the 
weak corrections and add them to the QED corrections [T7I . With this definition, both 
corrections are now AE independent. The main part of the QED corrections arises from 
these leading logarithms L e , originating from photons in beam direction. This leads to a 
large dependence on experimental cuts and detector specifications. We therefore use the 
structure function formalism ^H] and subtract the leading logarithmic 0(a) terms of the 
initial state radiation, a ISR,LL (s). After subtraction of these process- independent terms, 
only the non-universal QED corrections remain. This gives for the total cross section the 
final expression: 

(22) 



with 



a total (s) = 


a tree ( 


s) + a weak (s) + a QED (s), 


a weak (s) = 


a soft { 


» + = (<i-iO*f-fi.) 


a QED (s) = 


g.hard 


w-=(<i-wt-§*) 


a ISR ' LL (s] 




-L e f 1 dx ®(x)a tree (xs) , 

7T JO 






lim{5(l -ar)[- + 21og(e)] +9(1 



tree i 



a 



s) , (23) 



a tree (s) - a ISR > LL (s) , (24) 



(25) 

1 _|_ T 2 

x -e)-^?-}. (26) 
1 — X 

Further improvements would be to consider a more realistic electron spectrum and incor- 
porate beamstrahlung in the calculations. Due to their strong dependence on the actual 
experimental conditions, we do not include these effects. 



4 Numerical results 

For the numerical analysis, we concentrate on the production channels 

e + e" -> X1X2 and e + e" -> xlxl- 



S 



They are of special interest for future experiments because of their decay products and 
for kinematical reasons |2j. Due to the tree-level coupling structure, we study here two 
different scenarios: In the higgsino scenario the two lightest neutralinos are both nearly 
pure higgsinos and therefore the process is dominated by the s-channel Z° exchange. In 
the gaugino scenario with a bino and a wino as Xi an d X2 states, the selectron exchange 
diagrams play the most important role. In the following, we distinguish between the naive 
tree-level, the improved tree-level with the corrections to the neutralino masses m^o and 
the rotation matrix iV^ included, and the conventional weak and QED corrections to the 
improved tree-level as discussed in the last chapter. For the SM input parameters we use 
ot{m z ) = 1/127.922, m z = 91.1876 GeV, and m w = 80.423 GeV. 

4.1 Higgsino scenario 

For the definition of the higgsino scenario we use the following MSSM on-shell parameters 
in the convention |l(Jj : 

tan/3 = 10; // = -100 GeV; M 2 = 2M X = 400 GeV; M Q;L = M UtDfE = 350 GeV; A f = 400 
GeV; M^o = 700 GeV. This gives the one-loop corrected neutralino masses: 
Xi (94% higgsino): 87.8 GeV xl (97% higgsino): 110.0 GeV 
X° 3 (94% bino): 209.4 GeV xl (96% wino): 415.2 GeV 

In Fig. El we show the naive tree-level cross section for five different channels. The double 
higgsino production e + e~ — > XiXi with i — 1, 2 is highly suppressed due to the behaviour 
of the Z°x°iX j coupling. 

Numerical results for the radiative corrections to the X1X2 production are given in Fig. |3 
The total non-universal weak and QED corrections are in the range of -12% in the investi- 
gated parameter region and thus have to be taken into account in future experiments. In 
the case of X2X2 production, Fig. the small tree-level Z ^?,^ coupling leads to an en- 
hancement of the corresponding vertex corrections and to large box graph contributions. 
For the same reason, owing to the neutralino rotation matrix correction Nij the effect is 
also highly increased. Therefore, there is a big difference between the naive and improved 
tree-level cross section. 
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Figure 5: Neutralino pair production in the naive tree-level approximation with {full, 
dashed, dotted, dash-dotted, dash-dot-dotted} = {XiXi, X1X2, X2X2, X3X3, xlxl}- 

e+e~ Xi Xl e+e ~ Xi X2 




200 250 300 350 400 300 400 500 600 700 

y/a [GeV] y/s [GeV] 



Figure 6: Corrections to the X1X2 higgsino scenario. Left: The total cross-section in the 
naive tree-level approximation (dotted line), the improved tree-level (dashed line), and 
the full 0(a) corrected without ISR (solid line). Right: The full 0(a) without ISR (solid 
line), weak (dashed line) and non- universal QED (dash-dotted line) corrections relative 
to the improved tree-level. 

4.2 Gaugino scenario 

In the case of the gaugino scenario, we use as input the SPSla DR benchmark values 
[TS], defined at the scale Q = 454.7 GeV. With these values, we can calculate our on-shell 
parameters in a consistent way by subtraction of the corresponding counter terms, e.g. 
Mi = Mp*(Q) - 5Y n (Q), and obtain: 

tan/3 = 10.2; y. = 353.1 GeV; M x = 97.9 GeV; M 2 = 197.6 GeV; M A o = 393.6 GeV. 
In the sfermion sector, we only need the selectron mass parameters: = 198.0 GeV; 
Me = 138.0 GeV. For all other parameters, we can use the DR or on-shell values. The 
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Figure 7: Corrections to the X2X2 higgsino scenario. Left: The total cross section in the 
naive tree- level approximation (dotted line), the improved tree-level (dash-dotted line), 
with the weak corrections (dashed line), and the full 0(a) corrected one without ISR 
(solid line). Right: The full 0(a) without ISR (solid line), weak (dashed line) and non- 
universal QED (dash-dotted line) corrections relative to the improved tree-level. The 
dotted line shows the effect of the mass matrix corrections relative to the naive tree-level. 

differences are of higher order for our calculation. For the neutralino states we get: 
Xi (97% bino): 94.8 GeV xl (88% wino): 181.5 GeV 

Xl (99% higgsino): 360.3 GeV xl (88% higgsino): 377.4 GeV 

Note that the SPSla scenario is defined by DR parameter values. Thus the one-loop 
on-shell parameters given here can differ from those calculated in other renormalization 
schemes. The on-shell masses are of course the same up to higher orders. 
In Fig. |Hl we show the tree-level cross section for all three possible gaugino production 
channels and the higgsino X3X4 production. The double higgsino channel e + e~ — > 2x® 
with i=3,4 or mixed gaugino-higgsino channels are suppressed due to the given coupling 
structure. The full 0(a) radiative corrections for the X1X2 production, given in Fig 
are only in the few percent range because of the cancellation between the weak and 
QED corrections, especially near the threshold. While the QED corrections for the X2X2 
channel, see Fig. EH are also moderate and show a similar behavior to the previous case, 
the weak corrections strongly depend on yfs. For large y/s this can be studied in the 
so-called Sudakov approximation [20] . The corrections are -10% at ~ 750 GeV and even 
larger at higher energies. One reason is that the being mainly a wino, has also an 
11% higgsino component, which effects the weak corrections in a similar way to the X2 
pair production in the higgsino scenario. This results in a large negative correction for 
the sum of the QED and weak part. 
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Figure 9: Corrections to the X1X2 SPSla scenario. Left: The total cross-section in the 
naive tree- level approximation (dotted line) and with weak (dashed line), and full 0(a) 
(solid line) corrections without ISR. Right: The full 0(a) without ISR (solid line), weak 
(dashed line) and non-universal QED (dash-dotted line) corrections relative to the im- 
proved tree-level. 

5 Conclusions 



We have calculated the full one-loop elektroweak corrections to the neutralino pair produc- 
tion in e + e~ collisions. The chosen renormalization scheme can be used for the complete 
MSSM parameter space and all production channels e + e~ — > with i,j = 1,2,3,4. The 
process independent corrections to the neutralino mass matrix are included in the defini- 
tion of an improved tree level. We paid particular attention to an appropriate definition 
of weak and QED corrections. We extracted the non-universal QED corrections by sub- 
tracting the initial state radiation (ISR). The full one- loop corrections without ISR are in 
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y/s [GeV] y/s [GeV] 

Figure 10: Corrections to the X2X2 SPSla scenario. Left: The total cross-section in 
the naive tree- level approximation (dotted line) and with weak (dashed line), and full 
0(a) (solid line) corrections without ISR. Right: The full 0(a) without ISR (solid line), 
weak (dashed line) and non-universal QED (dash-dotted line) corrections relative to the 
improved tree-level. The dotted line shows the effect of the mass matrix corrections 
relative to the naive tree-level. 

the range of 5-20% and in some cases even larger and thus have to be taken into account 
in future linear collider experiments. 
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